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Abstract 

In previous publications ( J. Geom.Phys.38 (2001) 81-139 and refer- 
ences therein ) the partition function for 2-1-1 gravity was constructed for 
the fixed genus Riemann surface. With help of this function the dynam- 
ical transition from pseudo-Anosov to periodic (Seifert-fibered) regime 
was studied. In this paper the periodic regime is studied in some detail 
in order to recover major results of Kontsevich (Comm.Math.Phys. 147 
(1992) 1-23 ) inspired by earlier work of Witten on topological two dimen- 
sional quantum gravity. To achieve this goal some results from enumerative 
combinatorics have been used. The logical developments are extensively 
illustrated using geometrically convincing figures. This feature is helpful 
for development of some non traditional applications (mentioned through 
the entire text) of obtained results to fields other than theoretical particle 
physics. 
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1 Introduction 
1.1 Motivation 

It may not be exaggeration to say that dynamics of 2-1-1 gravity is just an in- 
terpretation of Nielsen -Thurston theory of surface homeomorphisms in terms 
of concepts known in physics [1-4]. Apparently, the reverse task of enrich- 
ing mathematics with concepts known from physics had been accomplished by 
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Wittcn [5]. His physical intuition had revolutionized combinatorial methods 
of algebraic geometry related to intersection theory on moduli space of curves 
[6-8]. In his famous paper [9] Kontsevich had provided needed mathematical 
justification to Wittcn's work. Since in both papers the initial object of study 
is 2 dimensional quantum gravity it is only natural to expect that the results of 
Witten and Kontsevich can be reobtained from more general 2+1 gravity model. 
The purpose of this paper is to demonstrate that this is indeed possible. 

It should be noted that the results of Witten and Kontsevich (W-K) rely 
heavily on the fact (proven by Kontsevich [9]) that the partition function of 2d 
gravity happens to be r function of the Kortcwcg-dc Vrics (KdV) hierarchy 
which is just a special case of more general Kadomtsev-Petviashvili (KP) hi- 
erarchy of nonlinear exactly integrable partial differential equations. Attempts 
to connect string theory with exactly integrable systems had been made ear- 
lier. Good summary of earlier efforts can be found in Ref.[10]. Recently, KP 
equations had been discussed in relation to dynamics of 2+1 gravity [11]. The 
content of Ref. [11] apparently has no connections with results of W-K. This 
does not exclude the possibility that such connection might exist and requires 
further study. In the present work we establish such connection (reduction) by 
rcobtaining W-K results from formalism developed earlier for 2+1 gravity [1-4]. 
This formalism is based on some properties of Riemann surfaces. From alge- 
braic geometry it is known that every Riemann surface S can be described in 
terms of the corresponding complex algebraic curve C [12] so that one can use 
interchangeably algebro-geometric and complex analytic language to discuss the 
same object, the Riemann surface. This,unfortunately, is not an easy task as 
was noticed by Looijcnga [7]. The existing difference in terminology is especially 
apparent when one is interested in the issue of compactification of the moduli 
space of Riemann surfaces (e.g. see Eq.(l.l) and related discussion) 

Moduli spaces Aig are connected with Teichmiiller spaces Tg of Riemann 
surfaces of genus (/ > 2 in such a way that Mg = Tg/T with F being the 
mapping class group of S. Both M.g and Tg have been discussed extensively 
in our earlier published papers on 2+1 gravity [1-4]. The concept of moduli 
had been initially proposed by Riemann who stated that isomorphism classes 
of closed Riemann surfaces of genus g > 2 are parametrized by 3g-3 complex 
parameters (or by 6g-6 real parameters). The space of allowed values of these 
parameters is effectively the moduli space Mg . More accurately, such moduli 
space is called coarse moduli space. It is complex analytic space of dimension 
3g-3 but not a complex manifold. Presence of singularities (to be discussed 
below) prevents this manifold from becoming a complex manifold [6,13]. One 
way of describing both Mg and Tg is through introduction of marking of <S 
(accordingly of C ). Marking of S can be made by introducing some n boundary 
components such that each of them is conformally equivalent to a punctured disk 
[14] . Accordingly, for algebraic curves we select n distinguished (smooth) points 
on the curve. Smooth means that the points are not located at the possible 
singularities of C. According to Deligne and Mumford [15] such singularities are 
regular double points which, in the case of traditional visualization of Riemann 
surfaces as closed surfaces with g holes, are associated with formation of nodes 
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Figure 1: Various types of degeneration of the Riemann surfaces 



as depicted in Fig.l. 

The compactification A4g_n of moduli space is acliicved by inclusion into 
Mg^n domains in parameter space associated with Riemann surfaces with nodes. 
In the language of algebraic geometry complex curves associated with such 
surfaces are called stable. Marked curves possess only finitely many authomor- 
phisms. It is known [6] that no stable curve can have more than 3g-3 nodes. This 
fact may be easily understood if we recall that each Riemann surface admits 
pants decomposition [2] so that 3g-3 nodes can be associated with degenerate 
geodesies whose lengths are squeezed to zero. For Riemann surface S instead of 
(or in addition to) marked points one can think about families of closed curves 
( lamination sets [4]) and their intersections. By properly defining intersec- 
tion number(s) one can measure the nontriviality of such sets [6]. In algebraic 
geometry such numbers are associated with cohomology classes known also as 
tautological classes ipi in H^{Mg^n, Q)- These numbers should not be confused 
with geometrical intersection numbers introduced by Thurston [16,17]. Roughly 
speaking such number i measures the number of intersections between a selected 
simple nontrivial closed curve a on <S and measured foliation i = a). This 
observation allows one to introduce equivalence relations between measured fo- 
liations. In particular, foliations Ti and T2 are projectively equivalent if there 
is number A G R+ such that i{J^i,a) = Xi{T2,ci). Following Thurston [17]and 
others [18] , we define the space of projective laminations VC{S) as space of 
equivalence classes of measured foliations. Such defined space forms a boundary 
of the Teichmiiller space 7^,n. The compactification of the Teichmiiller space 
T(iS) is defined now as 

'T{S)=T{S)\JVC{S). (1.1) 
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The set T(S) is closed ball of real dimension 6g—6+2n whose boundary is identi- 
fied with VC{S). Since the moduli space is the quotient of T(<S) by the mapping 
class group F , is clear that such defined compactification of the Teichmiiller 
space leads accordingly to the compactification of the moduli space. Such com- 
pactificatioon is apparently different from that of Deligne-Mumford. The differ- 
ence has some physical meaning. Indeed, in the traditional, algebro-geometric 
case, no dynamics is involved while in the Thurston's case the dynamics (that is 
"time evolution" ) is associated with dynamically formed 3-manifolds, e.g. see 
chapters 8 and 9 of Thurston's lecture notes [16]. These dynamical features are 
explained in physical terms and illustrated in our recently published paper on 
statistical mechanics of 2-|-l gravity [4] . 

Since in 2 dimensional quantum gravity there is no time evolution, naturally, 
there is no dynamical observablcs and theory is essentially topological. Witten 
[5] had suggested to study the topological averages of the type 

< T,, • • • rfc„ >g= [ < • • • ^^ (1.2) 

provided that 3g — 3 + n = ■ To understand the meaning of such an 

average it is instructive to consider the first nontrivial example: g = l,n = 1. 
In this case we obtain, 

< n >i= / V- (1-3) 

JMi,i 

Witten had demonstrated that < n >i= ^ • This result had been actually 
obtained earlier by Wolpert [19] who noticed that 

< Ti >i / LOw-p (1-4) 

where uj^-p is familiar Weil-Petersson Kahler (1,1) form which plays very im- 
portant role in both the theory of Teichmiiller spaces [20] and in string theory 
[21]. The dot above the equality sign has the following meaning. The integral 
Jj^ cOw-p had been independently calculated by Wolpert [22] and Penner [23] 

and is known to be equal to — so that the right hand side of Eq.(1.4) is equal to 

. The same result was obtained by Witten [5] who argues that it should be ad- 
ditionally divided by 2 because "the generic elliptic curve has two symmetries". 
Thus, evidently, computation of " tautological averages", Eq.(1.2), is effectively 
equivalent (up to numerical prefactors) to calculation of the Weil-Petersson 
volumes of M.g,n- Zograf [24] had calculated such volumes for n > 3 punctured 
spheres. Based on these results, Matone [25] had developed and solved non 
pcrtubativcly 2d Liouvillian gravity. Some of his results had been used and 
extended in the work by Kaufman, Manin and Zagicr [26]. These results had 
been further streamlined in recently published monograph by Manin [27]. 
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1.2 From Chern classes to Grassmannians 



To make an additional connection with physics, we would like to reanalyze 
just obtained results from the point of view of differential geometry of complex 
manifolds [21,28]. In this regard, a closed integral (1,1) form w on a complex 
manifold M determines the equivalence class of a line bundle A. Recall [21] that 
a line bundle A on manifold M is an assignment of one dimensional complex 
vector space Xz to each point z G M . Sections of Xz are some functions {z G 
Ba) assigning an clement of A^ to each point z. Vector spaces at different points 
should fit together and this leads naturally to the fiber bundle constructions 
made of transition functions (pa^ between different coordinate charts Ba- A 

I 1 2 

metric on Xz is a set of positive functions Qa such that ga = |<?^a^| 5/3- The 
covariant derivative is defined now as 

The bundle A is holomorphic if in addition 

V,/„ = . (1.6) 

With help of such defined covariant derivative the curvature tensor F of A is 
obtained in a standard way [28] with the result: 

F=^ln,.. (1.7) 
Finally, it can be shown [21,28], that 

to = -^Fdz A dz (1.8) 

and that, actually, such defined (1,1) form coincides with the first Chern class 
form ci(A). In our case M = M.g,n and the line bundle is actually a cotangent 
bundle T*{Mg^n) made of quadratic differentials [Wolpert, positive] .Wolpert 
[29,30] had demonstrated that 

\iOn,-p = ^F = Cx{X). (1.9) 

Substituting this result into Eq.(1.4) we obtain, 

< Ti >i= / cip(A), (1.10) 

JMi,i 

where cip(A)=Ci(A)/2 . In such notations this result coincides with Eq.(1.7) of 
Witten's paper [5] (for d= 1). The reason of introducing the factor of 1/2 can 
be explained as follows. According to Wells, Ref.[28], Ch-r 6, the traditional 
construction of Hodge line bundles for manifolds requires that 

= integer = c. (l-H) 
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Moreover, if c is an integer, it is actually equal to the Euler characteristic of 
the manifold. In our case, Mi^i is an orbifold and for the orbifolds the Euler 
number can be rational number according to Thurston [16] . Wells explains how 
to construct the Hodge bundle for any Kahlcrian manifold/orbifold. To this 
purpose, if the constant c in Eq.(l.ll) is not an integer, it is sufficient to rescale 
the form u; : O = c~^ci; , so that the rescaled form is the Hodge form. Since the 
moduli space A4g.n is an orbifold, Eq.(l.lO) is acceptable and, hence, with such 
remarks, coincides with Eq.(1.7) of Witten's paper (for d = 1), Ref.[5]. 

Obtained results possess additional very important physical information 
allowing us to make direct and simple connections with Grassmannians and, 
hence, with exactly integrable systems, using arguments different from that of 
Witten [5] and Kontsevich [9]. 

Let G'~^{d,n) be complex Grassmann manifold. It is known that 

where U{n) belongs to the unitary Lie group of n x n matrices. Following 
Ref. [31,32] consider now the classifying space BU{n) defined by 

BU{d) = lim (1.13) 

ra^oo U [d)U (n) 

The name classifying comes from the following observation. Suppose there is a 
map / from a base space M to BU{d), i.e. f: M ^ BU{d). Construction of 
such a map for the case of orbifolds has been developed by Baily [33] and, in 
particular, for M.g,n by Wolpert [29]. The n- vector complex bundle A over the 
base space M can be expressed as a pullback of the standard vector bundle ^ 
over BU{d) . That is A = f *{£,). The i-th Chcrn class c* of the vector bimdle 
A can then be calculated simply by the pullback f*{ci) G H^^{M,R) of the 
element Cj in the 2i-th cohomology group of BU (d). G^{d, n) can be embedded 
in the complex projective space. This can be achieved for any n and the limiting 
case n ^ 00 is known as Sato Grassmannian in the theory of KP equations [34] . 
We would like now to describe such an embedding in some detail. 

Recall that the complex projective space P" is defined as P" := (C"+^ — 
{0})/ ~ with equivalence relation ~ defined as follows. If a point P eP" is given 
by an n + 1 tuple (z(0), z(n)), then another (n + 1) tuple {z/{0), z/{n)) 
defines the same point P gP" if there is a nonzero number c such that 
z{i) = cz/{i), yi , i = — n.The system of standard coordinates {Ui enables 
us to define a manifold structure on P" : 



and 



[/, :={(z(0),...,z(n))|z(i)9^0},* = l 



.TT ^C-n ._ fZ{0) 2:(i-l) z(i+l) z{n) 
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With help of such map Ui can be identified with the afiine space C" and 
the space P" can be made from different coordinate patches Ui with help of 
transition functions (pa/, as discussed before. A linear space jC in P" is defined 
as the set of points P={z{0), z{n)) of P" whose coordinates satisfy a system 
of linear equations 



r. 



= (115) 



a = 1, (n — (i). The space C is d-dimensional if [n — d) x (n + 1) matrix of 
coefficients \baj\ has a nonzero {n — d) x (n — d) minor. In this case there are 
d+1 points Pi = {zi{0), ...,Zi{n)) in C (i = — rf) which span C . Naturally, 
£ is a line ii d = 1, a plane ii d — 2 and a hypcrplanc if d > 2. Wc shall call 
these planes as d-planes following [35] . d-planes in P" can be represented by 
the points in the projective space P^ whose dimension N is given by 

''= idlmn-d)r '- ^'-''^ 

To this purpose, let us fix a rf-planc in P" and pick d+l points Pi =(zi(0), Zi{n)) 
which span £. Using these points let us form (d + 1) x (n + 1) matrix \pi{j)] 
with < i < d and 0< j < n. Let jo, ■■■■,3d be a sequence of integers with 0< 
3(i ^ n and let p(joj ■■■ijd) denote the determinant of (rf + 1) x (rf + 1) matrix 
IPiiJf})] ■^^itih i,P = 0, There will be A^+l determinants of such type and 
at least one of them is nonzero by requirements of linear algebra. Hence, in 
view of Eq.(1.14), wc conclude that we can use these determinants to determine 
a point in the complex projective space P^. The coordinates of this point are 
called Pliicker coordinates of C in P^ and such an embedding of the com- 
plex Grassmannian manifold (of c?-planes in P" space) into complex projective 
space P^ is called Pliicker embedding. Not every point in P^ arises from em- 
plane in P". Pliicker coordinates p{jo, ■■■,jd) obey the following set of (Pliicker) 
equations 

Ed+l 
. {-iyp{3o,.:,jd-ikj)p{kQ,...,kj,...,kd+i) = Q. (1.17) 
3=0 

Here jo,...,jd,-\ and fco, ., fc^+iare sequences of integers with 0< jp , k^ < n 
with kj meaning that the integer kj has been removed from the sequence. 

As it is shown by Miwa et al [34] Pliicker coordinates represent the location 
of tau function of the KP hierarchy inside the Grassmannian while Pliicker 
equations are in one to one correspondence with the Hirota bilinear equations. 
Hence, the connection between the averages given by Eq.(1.2) and KP (or, more 
exactly, KdV) hierarchy naturally follows. The determinants, which are points 
in P^, have probabilistic meaning which is discussed below. To this purpose 
we have to introduce some additional concepts. 

1.3 From Grassmannians to Schubert varieties 

Consider a sequence of subspaces (cellular decomposition), i.e. Aq C Ai C 
... C of a fixed space, e.g. P", each properly contained in the next, whose 
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dimensions dim Ai = ai , provided that < aq < ai < ... < < n . Such 
construction is called flag. Let now ^{Aq, A^) be the subset £ of G^{d,n) 
consisting of all d-planes satisfying dim (jCdAi) > i for i=0, ...,d. Thus, if dim 
Ai = i Vi , then il{AQ, ...jAd) is made of a single d-plane, while if dim^dj = 
n-d + iyi , then n{Ao, A^) = G^{d, n). 

Definition 1.1. Q,{Aq, Aj) is called Schubert variety corresponding to 
the flag j4o C Ai C ... C Ad . VI{Aq , . . . , Ad) defines a homology (actually, 
cellular homology [32] ) class in the homology ring H^^{G'^{d,n);Z). 

Definition 1.2. Homology class in n); Z) is called Schubert cy- 

c/e.Becausc the homology class depends only on the integers =dm\Ai , it is ap- 
propriate to write ^{Aq, ...,Ad) = ri(ao, ■■■ad) where 0< ao < oi < ... < ad < 

It can be shown [12] that the product of any two Schubert cycles can be 
uniquely expressed as a linear combination of other Schubert cycles. This obser- 
vation is central for development of Schubert calculus [12,35,36]. 

Remark 1.3. Gepner [37] had demonstrated that all results of rational 
conformal field theories can be actually obtained from physically reformulated 
Schubert calculus. Additional physical refinements of these ideas can be found 
in the work by Witten [33] . 

H^{G'^{d, n); Z) is generated by the special Schubert cycles given by 



for i = 0,1,.. .,n — d. These results allow to prove [12,35,36] the following 
theorem of central importance for the whole development presented in the rest 
of this paper. 

Theorem 1.4. For all sequences of integers < ao < ai < ... < ad < n 
(which we denote as a ) the following determinantal (Giambelli's-like) formula 
holds in the homology ring H^{G'^{d, n); Z) : 



with \a{ai + j — i)\ being a determinant made of special cycles. 

Remark 1.5. The name Giambelli comes from the fact that structurally 
the same expression exist for the Schur polynomials Sa which was discovered by 
Giambelli [39-42]. The Schur polynomials are characters of the general linear 
group on symmetrized complex linear vector space E- [39-42]. In the light of 
results presented earlier we can associate the determinant p{jo, ■■■^jd) with that 
for special Schubert cycles, Eq.(1.19), so that the Schur polynomials are tau 
functions of KP hierarchy. The formal correspondence between the Schur poly- 
nomials and Eq.(1.19) is not coincidental. It can be proven [42] that, actually, 
there is an isomorphism © between Sa and f2(a). 

Remark 1.6. The topological meaning of Eq.(1.19) had been clariflcd by 
Portcous [43] and also by Horrocks [44] and, later, by Carrel [45]. All these 
results can be actually deduced directly from much earlier fundamental papers 
by Chern [46] and Ehresmann [47]. 



a{i) = fl{i, n— d+1, ...,n) 



(1.18) 



V,{&) = \a{ai + j — i)\ ,0<i,j<n — d 



(1.19) 
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We would like to describe briefly these results since they are essential for 
correct physical understanding of the meaning of intersection numbers. To 
this purpose let us observe that for fc-cycle A and n — k cycle B on an n- 
dimensional manifold M the Poincare ' duals of these cycles arc closed n — k 
and k differential forms ip and ip respectively so that the intersection number 
*{A- B) of these cycles in homology is equal to the wedge product of these two 
forms in cohomology [12], i.e. 

*{A-B) = y <^A V- (1.20) 

M 

This definition can be extended to describe the intersection of subvarieties V and 
W oi a, complex manifold M and it is possible to prove that Schubert calculus 
is just a special case of such more general algebra known in the literature as 
Chow algebra [8] (also as Chow ring [41]). Next, we need a notion of a divisor. 
If the manifold M can be decomposed as 

M = ViU- ■ - UVm (1.21) 

then one introduces 

Definition 1.7. Divisor on M is locally finite formal linear combination 

D = ^a,T/, (1.22) 

of irreducible analytic hypersurfaces 1^ of M. 

Definitions of "local finitness" and " irreducibility" are given in Ref.[12] on 
page 130. The constants are normally some integers. For G'^{d,n) Schu- 
bert cycles provide the desired decomposition of the Grassmannian [12]. The 
following theorem is of central importance [12]. 

Theorem 1.8. The Chern class ci(A), given by Eq.(1.9), of the line bundle 
A represents the Poincare' dual of the fundamental homology cycle carried by 
the divisor D, i.e. 

^yFAV = ^a.y^ (1.23) 

M Vi 

for every real closed 2n — 2 form t/j. 

Consider now some implications of this theorem. First, in the case if M is 
compact Riemann surface, a divisor on M is just a finite sum 

D = J2niPi (1.24) 

of points Pi G M with multiplicities n^. Combining Eqs.(1.23) and (1.24) 
produces the Poincare'-Hopf index theorem 

^jF = Y,ni=x{M). (1.25) 

M * 
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This theorem played central role in our earlier works of 2+1 gravity [1,2]. This 
observation clarifies both the meaning of numbers rij and points pi : for vector 
(or line) fields these points are associated with singularities of the field. These 
singularities had been interpreted as masses. 

Remark 1.9. Generalization of this result to higher dimensions had been 
developed by Chern and Weil (CW). The up to date exposition and general- 
ization of their results can be found in Ref.[48] and, in principle, provides an 
opportunity to describe 3+1 gravity in a way similar to that developed for 
2+1 dimensional case. Such approach to gravity is very close in spirit to that 
originally suggested by Regge [49]. 

We need this observation in the present context as well for the following 
reasons. By definition, a subset V of an open set U G C" is an analytic variety 
if for any p E U there exists a neighborhood U of p in [/ such that V O U' is 
common zero locus of a finite collection of holomorphic functions {fi} on U . 
In particular, V is called an analytic hypersurface if locally V is zero locus of a 
single nonzero holomorphic function /, i.c.V = {f{z) = 0} in the neighborhood 
of e Accordingly, for decomposition, Eq.(1.21), Vi = {fi{z) = 0} with Vi 
irreducible at . All these facts lead to the following 

Definition 1.10. An algebraic variety V C P" is the image in P" of zero 
locus of a collection of homogenous polynomials defined in C"+^, i.e. V={ 
J^i{zo,...,Zn) = 0). 

For the line bundle A on M it is possible to associate sections with if 
M can be embedded into P"(and, in our case, it can be embedded since we 
had mentioned already that A = /*(^) with ^ e BU{n)). Suppose now that for 
some point z* G M of n-dimcnsional complex manifold m sections arc linearly 
dependent. Then, according to CW theory [48,50], this fact can be written as 

/ai A---A/,^ =0. (1.26) 

This equation is multidimensional analogue of the Poincare' condition (e.g. read 
Remark 4.2. of Ref [2]) for the singularities of the line/ vector fields on surfaces. 
The degeneracy set is set of all points z* for which condition given by Eq.(1.26) 
holds. In the most general case it is m — 1 dimensional submanifold D of M. 
Consider now a cycle a on M of dimension r lesser than that of D. Let such a 
cycle meet D transversely at the point z* eD then, f{a) will meet some Schubert 
cycle r2(a) also transversely at the point f{z*) of the Grassmannian so that the 
intersection number of /(a) with f2(a) at f{zl) will be the same as that for a 
meeting D on M. Taking into account Eq.s(1.20) and (1.23) and also Theorem 
1.8., we conclude that Cr(A)(a) = *(a,D). With little additional work it can be 
shown [12] that 

n{a) = \ca,+j-i{\)\ ,0<i,j<n-d, (1.27) 

to be compared with Eq.(1.19). This result was obtained by Porteous [43] and 
is known in the literature as Porteous formula [12,41]. 
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1.4 From Schubert veirieties to directed random walks 



Porteous formula can be seen as a special case of much more comprehensive 
result of Weil (and developed by Chern ) known as Weil homomorphism. In 
view of this homomorphism, the result given by Eq.(1.27) can be reinterpreted 
as Schur polynomial (e.g. read discussion following Eq.(1.19)) and, hence, as r 
function of KP hierarchy. Schur polynomials Sa originate from the known [40] 
identity for indeterminates {xi} 

(Xi + • • • + XmT = /^Sa(xi, Xm), (1-28) 
X\-n 

where use of the notation A h n is meant to say that A is partition of n. Partitions 
are best represented by the Young tableaux. Accordingly, the factor denotes 
a number of standard tableaux of a given shape A . 

Eq.(1.28) can be rewritten in a slightly different form as follows 

(a;i + ...+xX= E f^^-'-C-- (1-29) 

, ^ ^ mi!m2!...m„! 

provided that M = mi + • • • + m^. This observation allows us to write Schur 
polynomial Sa in the form 

SA(a;i,...,a;„)= ^ ^A.m^^ • • • C (1-30) 

m=(mi,...,mn) 

with K\^rn being some coefficients (Kostka numbers) [42]. 

Eq.(1.30) can be given probabilistic meaning in terms of the directed random 
walks. Indeed, following Ref.[42,51], consider planar lattice. On this lattice 
consider a directed path P from (a,l) to (b,N). The information about this path 
can be encoded into multiset HoTyiV) of y— coordinates of the horizontal steps 
of P. Define 

w(p)= n (1-31) 

i=Hory{P) 

To facilitate reader's understanding, we illustrate these ideas on Fig. 2. In this 
figure iIory(P)={l,3,3,6} and, accordingly, w(P) I. Next, we need to 

extend this result to an assembly of directed random walks ("vicious" random 
walkers in terminology of Fisher [52]). That is we need to consider the products 
of the type w(Pi) • ••w(P/s) =W(P). Finally, the generating function for an 
assembly of such vicious walkers is given by 

hb-a{Xu-,XN) = Y^(P) (1.32) 
P 

where W(P) is made of monomials of the type a;™^a;^^...x™" provided that 
mi + ■ ■ ■ + niN = b — a. The following theorem is proven in Ref.[42,51] 
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Figure 2: A typical directed random walk 



Theorem 1.11. Given integers < ai < ... < Ofc and < bi... < bk , let 
j be the k x k matrix 



then, 



Mij = hb^-ai{xi,...,XN), 



detM=^W(P), 



(1.33) 



(1.34) 



where the sum is taken over all sequences (Pi,...,Pfe) =P of non-intersecting 
lattice paths Pi : (a,, 1) {bi, TV). 

Corolleiry 1.12. Put now = i and bj = Xi+j in Eq.(1.33) provided that 
^ < i,j < k with A being partition of N with k parts then, detM=SA(a;i, ...,xjv). 

Theorem 1.11. provides desired connection between the vicious random 
walkers and the r function of KP hierarchy. 



1.5 Organization of the rest of this paper 

In Section 2 we provide some facts from the theory of random vicious walkers 
using results of Fisher [52], Huse and Fisher [53] and Forrester [54] .We argue 
that these results can be obtained also with help of the Bethe ansatz method 
applied to one dimensional non ideal Bose gas. Such observation is helpful for 
developing connections between the Yang-Baxter equation, symmetric functions 
and Schubert polynomials. The obtained Bethe ansatz wave function is reinter- 
preted in terms of the Gaussian unitary ensemble of random matrices. Using 
some results for random matrices summarized by Mehta [55] and more recent 
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results by Tracy and Widom [56] (and Forrester [57]) wc discuss relations be- 
tween the Kontsevich matrix Airy integral and that coming from the Gaussian 
unitary ensemble. We argue that both integrals are tau functions of KP hierar- 
chy. Moreover, wc demonstrate that these integrals arc actually tau functions 
for KdV hierarchy of equations (that is they are only a special case of KP hi- 
erarchy) and, hence, both can be used as solutions of the W-K model. The 
arguments of Section 2 are too general. They do not contain explicit reference 
to 2+1 gravity, moduli space, etc. This deficiency is corrected in Sections 3 
and 4 which contain new combinatorial proof of the main identity obtained by 
Kontsevich [9]: 

^ {2d, - 1)!! _ 2-^(r)2g(n 1 

2^ <Td,---Td^>g \_l 2rf^^i - 2^ \AutiG)\ 11 \(e)+Xi(e)' 

Sdi=3a-3+n i=l ^^ GeTg^n ' ^ ^' ee£(G) ' ^^^> 

(1.35) 

This identity connects the observables of topological quantum gravity, Eq. (1.2), 
with averages of the random matrix (Kontsevich) model associated effectively 
with the ribbon graphs Tg^n representing the combinatorial moduli space A^g"™*" 
of marked Ricmann surfaces of genus g. To avoid repeats of technical details 
describing random matrices, ribbon graphs, etc., discussed in Sections 3 and 4, 
we only notice that V(r) in Eq.(1.35) represents the total number of vertices 
and £{r) the total number of edges of the ribbon graph Tg,n- Mathematically, 
the identity, Eq.(1.35), is just the statement that there are different but equiva- 
lent ways to present partitions discussed in Section 1.4. In section 3 we connect 
these partitions with some results coming from Nielsen- Thurston theory of sur- 
face automorphisms. As it was demonstrated in Refs.[l-4], this theory provides 
natural mathematical framework for description of dynamics of 2-|-l gravity. In 
these references the canonical (fixed genus g) partition function for 2+1 gravity 
was obtained. This partition function was used for description of the dynamical 
transition from the pseudo-Anosov to Seifert-fibered (periodic) regime (phase) 
of 2+1 gravity. We remind our readers about these concepts within the bound- 
aries of Nielsen- Thurston theory. In Sections 3 and 4 we argue that the W-K 
partition function is relevant for the Seifert-fibered phase. This phase was not 
discussed in Refs.[l-4]. We should warn our readers that complete description of 
the Sefert-fibered phase requires more than just the W-K model. For instance, 
Kulkarni and Raymond [58] had found a very interesting connection between 
the Seifert-fibered and anti-de Sitter 3 manifolds which are just Lorentz mani- 
folds of constant negative curvature. Such anti-de Sitter 3 manifolds had been 
recently classified by Francois [59] . In spite of the fact that they have received 
considerable attention in physics literature [3] recently, the complete description 
of this phase ( which by the way contains all known crystallographic groups and 
much more [60] ) in physical terms remains challenging research problem. 

Thus, the results of Refs.[l,4] for the partition function of 2+1 gravity should 
be extended to reach an accord with results of W-K model. To this purpose, fol- 
lowing Nielsen [61], dynamics of the Riemann surface homeomorphisms should 
be lifted to the universal cover, e.g. to the Poincare' disc model of H^. In such 
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model the set of geodesies (geodesie lamination) is represented by the set of 
non intersecting arcs (i.e. circular segments whose ends lie on 5'43).The com- 
binatorial arrangement on the disc is described by the Catalan numbers C„. 
These numbers had been used earlier in our work, Ref.[l] for construction of 
the partition function. In this paper we use this combinatorial data in order to 
establish several important bijections:!) from the set of arcs to the set of Dyck 
paths, Fig. 5; 2) from the set of Dyck paths to the set of parallelogram polyomi- 
noes, Fig. 6; 3) from the set of polyominoes to the set of Young tableaux and, 
finally, 4) from the set of Young tableaux to the set of vicious walkers, Figs. 7 
and 8. This chain of bijections is needed in order to bring our earlier obtained 
results for the partition function in accord with the l.h.s. of the Kontsevich 
identity, Eq. (1.35). To bring this partition function in accord with the r.h.s. of 
this identity the ribbon graphs need to be constructed. This is discussed in 
Section 4. Unlike Kontsevich [9] and others [62-64 ], we do not use quadratic 
differentials (discussed earlier in our papers, Refs.[l,2] ) for construction of these 
graphs. Our method is based on use of known various equivalent geometrical 
ways to describe the combinatorics of Catalan numbers [42] . The edges of thus 
constructed ribbon graphs are replaced by the paths of vicious walkers so that 
in the end using combinatorics of the Young tableaux the partition function for 
assembly of vicious walkers acquires the form of the r.h.s. of Eq.(1.35). 

Section 5 is included in this work for several reasons. First, it discusses the 
issue of universality of the W-K model from the point of view of dynamical 
systems theory. By universality we mean the fact that partition functions of W- 
K and many related models to be discussed briefly in this section are solutions 
of the KdV hierarchy. The reasons for this universality can be traced back to 
the very basic properties of quasiconformal transformations known already to 
Ahlfors long time ago [65]. Although in Ref. [3] these transformations had been 
discussed extensively, in this paper more up to date information is presented to 
make it relevant to the results obtained in earlier sections. This includes some 
facts about the Thompson and the Ptolemy groups, about their connections with 
binary trees (and, hence, with Catalan numbers) and their relations with the 
universal TeichmiiUer and moduli spaces, etc. Second, it discusses briefly some 
potential physical and biological applications of the obtained results. Finally, 
it discusses connections with Frobenius manifolds, self-dual Einstein equations, 
etc., thus leaving many problems open for further study. 

2 From vicious walkers to Kontsevich model via 
Gaussian unitary ensemble of random matri- 
ces 

Following Forrester [54] and Fisher [52], we would like now to formulate a lock 
step model of vicious walkers. Incidentally, this model is just many walks 

generalization of the directed polymer model considered in our earlier work, 
Ref. [66] .The continuum limit of the distribution function for this directed poly- 
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mer produces the Euclidean version of the Dirac propagator for particle whose 
mass is associated with bending probabilities to be discussed below. In case of 
many walkers the fermionic nature of the Dirac "particles" (walks) imposes a 
sort of Pauli principle which forbids two walks to intersect. This is characteris- 
tic to all quantum many body problems where all "particles" live in the same 
"world time" [67]. In the theory of Brownian motion each walker is allowed to 
have its own world time [68] so that, accordingly, one can have many world 
times quantum mechanics. We shall refrain from discussion of these options re- 
ferring interested reader to current literature [69]. The lock step model assumes 
just one "world time". Technically this means the following. 

We consider a square lattice where x coordinate is assigned for "space" 
while the y coordinate is assigned for "time" . If p walkers are labelled in linear 
sequence along x axis so that one has 

Xi < X2 ■ ■■ < Xp (2.1) 

at "time" t = 0, the same inequalities should hold for all subsequent times. The 

walkers start either on the even or on the odd numbered sites on the x axis. At 
each tick of the clock each walker moves either to the right or to the left (along 
lattice diagonals) with equal (bending) probability w. The probabilities, in 
general, may not be equal and are associated with masses of the Dirac particles 
as discussed in our work, Ref.[66]. Because of imposed initial condition, no 
two walkers can occupy the same lattice space at any time. This circumstance 
makes walkers "vicious". Let Xq = (a;i^o,...,2^p,o) be the initial configuration of 
such vicious walkers and x=(xi, ...,Xp) be the final configuration at time t. To 
calculate the total number of walks starting at f = at xq and ending at ta.t 
X we need to know the probability distribution Wp (xq x; t) that the walkers 
proceed without passing, i.e. maintaining the inequalities 

Xj-i{t') < Xj{t') , j = 2,3, ...,p iovO<t' <t (2.2) 

from their initial configuration at time t = to their final configuration at time 
t. It is of interest to study this problem for large times: t ^ oo . In this case 
it is possible to make a transition to the continuum limit in our study of the 
probability distribution. In doing so, we shall follow the arguments of Huse 
and Fisher [53]. To begin, we need to recall that in the continuum limit the 
probability distribution Wp(xo — > x;t) for p independent random walkers is 
known to be 

W^i^o - x;i) = exp{-^i^^}/(27rl?t)^'/2. (2.3) 

The diffusion constant D sets up the scale since, as usual, <(xj — .Xj.o)^ >= Dt. 
Because of use of a single world time for all walkers the above distribution 
function can be also viewed as distribution function for a single walker in p- 
dimcnsional Euclidean space. The restrictions given by Eq.(2.2) impose addi- 
tional constraints that such a walker must not cross any (hyper) planes described 
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*U,-x,=o+=2c*,l<i<i<p. (2.6) 



by the set of equations xi = X2,X2 = .T3, ...,Xp^i = Xp. In view of the rcsuhs 
of Section 1, e.g. see Fig. 2 and discussion related to it, the assembly of these 
planes forms Grassmann manifold. This statement can be clarified further by 
considering the following example. Following Gaudin [70] consider Schrodinger 
equation for one dimensional Bose gas of p particles interacting via point-like 
pairvise interaction potential. The dimensionless form of the Schrodinger equa- 
tion for such particle system is given by 

- E U + 2c E ^(^^ - ^^O* = (2-4) 

This equation is equivalent to the boundary value problem of obtaining the wave 
function of the equation 

-Ap* = (2.5) 

where Ap is just p-dimensional "free" Laplacian, and the wave function ^ is 
subjected to the set of constraints: 

_a d_ 

^ dxi dxj ^ 

The condition ^'=0 on the hypcrplanes is achieved in the limit c=oo according 
to Gaudin [70]. It can be shown, that for any c solution of the corresponding 
quantum mechanical problem is obtained with help of the Bethe ansatz method 
[70]. Moreover, any problem solvable by the Bethe ansatz method is essen- 
tially of the type just described as had been demonstrated rigorously by Gutkin 
[71]. Therefore, not surprisingly, that there are deep connections between the 
classical exactly integrable systems of KP type and the quantum mechanical 
exactly integrable systems solved by the Bethe ansatz method [72]. The Hecke 
algebra leading to the Yang-Baxter equations providing mathematical justifica- 
tion of the Bethe ansatz method is coming from some particular representation 
of the symmetric group [73] and, hence, is connected with Schur and related 
polynomials. Some additional details can be found in Ref. [74]. 

The mathematical problem posed by Eq.s(2.5) and (2.6) can be equivalently 
formulated as problem about properties of random walk inside p— dimensional 
caleidoscope. That is we are looking for solution of an eigenvalue problem for 
"free " Laplacian compatible with some reflection group[70,71]. In the simplest 
cases these are just subgroups of the symmetric group Sp made out of even and 
odd permutations. In view of this, we obtain, 

Wpi^o ^ x;t) = ^ e(5)W^p(5Xo ^ x;i), (2.7) 
geSp 

where £(17) = ±1 depending upon the symmetry of permutation (even or odd). 
Taking Eq.(2.3) into account, this result can be rewritten equivalently as 

Wpi^o - x;t) = ^p(xo,x;t) ^"P^-f (2.8) 
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where 

C/p(xo, x; t)=J2 <9) exp [{x-gxo)/Dt] . (2.9) 

Some short calculation explained in Ref. [53] produces 

C/p(xo,x;i) ~ consiA(x)A(xo)/ (Di)"*" (2.10) 

with const = l/l!2!...(p— 1)! , Up = ^vip— 1) and A(x) being the Vandermonde 
determinant: 

A(x) = n(a;i-ar,). (2.11) 

Remark 2.1. From standard texts in probability theory, e.g. see Ref. [75], 
it is known that non normalized expression for the probabiUty (xq x; t) 
is the long time limit of the formula providing the total number of walks of 
n steps (since t ^ n) from point Xq to point x. Accordingly, Eq.s (2.7)- 
(2.10) provide the total number of nonintersecting directed walks and, hence, 
Wp(xo ^ x;t) w det M as shown in Eq.(1.34). 

It is convenient to assume now that x=xo. Then, upon rescaling, the follow- 
ing result holds: 

Wp(x = xo;t) = Pp{xi, ...,Xp) = const x exp(-x^)A^(x). (2-12) 

This is just the probability distribution of eigenvalues of random matrices from 
the Gaussian unitary ensembe [55]. Following Dyson [76] define the n— point 
correlation function by 

Rn{xi,...,Xn) = J ■ ■ ■ J Pp{xi,...,Xp)dXn+l ■ ■ ■ dXp. (2.13) 

— oo — oo 

Using this definition the connected n-point correlation functions T„(a;i, a;„) 
are defined in a usual fashion, e.g. 

Ti{x)=Riix) 

T2{xi,X2) = -R2{X1,X2) + Rl{xi)R2{X2), 

etc. Using method of orthogonal polynomials [55] it can be demonstrated that 

Rn = det[Kp{xi, Xj)]iij=i^,„^ri (2.14) 
where the kernel Kp{xi,Xj) is given by 

p-i 

Kp{x,y) = ^ipk{x)ipk{y) (2.15) 
fe=o 
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with functions (fk{x) depending on the random matrix ensemble used. In the 
case of Gaussian unitary ensemble they are given by 

Mx) = (2'=fc!\/^)-5 exp(a;V2)(-^)^' exp{-x^). (2.16) 

The connected n-point correlation function Tn{xi, ...,Xn) can be neatly repre- 
sented using the kernel Kp{x,y) as follows [55], page 92,: 

Tn{xi, ...,Xn) = ^Kp{xi,X2)Kp{x2,X3) ■ ■ ■ Kp{ ), (2.17) 

p 

where the sum is over all (n — 1)! distinct cycUc permutations of indices (1, n). 
For p — » 00 the kernel Kp(x,y) can be calculated and in terms of the rescaled 
variables it was obtained independently by Tracy and Widom [56] and by For- 
rester [57]: 

K^X,Y) = ^&^^^^=4^n^, (2.18) 
where Ai{y) denotes the Airy function 

CO 

Ai{y) = j exp(i(a;^/3 - xy))dx (2.19) 

— oo 

and the prime denotes differentiation with respect to its argument. Following 
Kontsevich [9] we define now the matrix Airy function analogous to the "scalar" 
case 

A{Y) = J exp(i tr(XV3 - XY))dX (2.20) 

where X and Y are hermitian N x N matrices for some N. After some compu- 
tation Kontsevich obtains, 

^(Y) = ,2.)"<"-"^i5^j4^, (2.21) 

det(l^-' ) 

where 

A^i-^\yi)= j dxxi-^eMii.x^l^-xyi))=(^-^^ Ai{yi). (2.22) 

— OO 

For 2, = 1 — 2 we obtain, 

A{X,Y)/{2t:YI^ = K{X,Y), (2.23) 

that is the Tracy- Widom kernel and the Kontsevich Airy matrix integral are 
practically identical. Naturally, it is of interest to find out if this result will 
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hold for i, j > 2. Comparing Eq.s. (2. 14), (2. 17) and (2.23) we conclude, that 
Eq.(2.17) should be considered as a likely candidate for further treatment. This 
conclusion is in accord with recent results of Okounkov [77]. We shall use some 
of his results below in Section 3 while in this section we would like to discuss 
different approach. To this purpose, following Mehta [55] , let us notice that 
correlation function i?„ , Eq.(2.13), can be presented in the following form 

m 

= ^(-1)"— H Kp{Xa,Xb)Kj,{Xb, Xe) • • • Kp{Xd, Xa) (2.24) 

p 1 

where the permutation P is a product of m exclusive cycles of lengths hi,h2, ■■■,h„ 
of the form {a^b^c^---^d^a), ^™ hj = n. Comparison between 
Eqs.(2.17) and (2.24) indicates that if, say, the connected n— point correlation 
function T„ is tau function of KP hierarchy, then correlation function i?„ should 
possess this property as well. 

To prove that both i?„ and T„ are indeed tau functions several steps are 
required. First, we would like to reconsider Eq.(2.21) in the light of subsequent 
refinements of Kontsevich work in physics literature. Following Di Francesco 
[78] the Kontsevich integral 6 n (A) is given by 



e;v(A) = 



(2.25) 



where A is diagonal N x N real matrix with elements A = (Ai, A^v) along the 
diagonal and M being N x N Hermitian matrix. After some calculations this 
integral can be brought to the following form: 



ejv(A) 



7V-1, 



\z,Dz,D^z,...,D 
\1, A, A^, A-'^"-^] 



(2.26) 



where 



^(A) 



dm\ — exp 



*<-6- + — 



3 ' 



(2.27) 



and 



D = X + 



1 



1 d 



(2.28) 



2A2 XdX' 

The Vandermonde determinant A(A), e.g. see Eq.(2.11), is written in the present 

case as 



A(A) = |1,A,A2,...,A 



N-l I 



(2.29) 



so that the expression in the numerator of Eq.(2.26) is also a determinant. 
Second, since z{X) is a solution of Airy's equation 



(D^ - X') z{X) = 



(2.30) 



19 



written in a somewhat unconventional form, it is possible to replace terms of 
the type D^Pz{X) in the determinant of Eq.(2.26) by A^Pz(A) and, analogously, 
£>2p+i^(A) by A2p+i0(A) where 

^(A) = \dz{X). (2.31) 
This allows us to rewrite ©jv(A) in the following form 

|„iV-l„ ^iV-2= I 

with X = 1/X. Third, if the asymptotic expansions of z and z given by 

0(A)=^CfeA-3'= (2.33) 

fe>0 

z{X) = J2dkX-^'' (2.34) 

fc>0 

with known coefHcients Ck and dk are substituted into Eq.(2.32) it acquires the 
following final form: 

e.(A)= E n«r"^^^ — J-^x--^ 11 • ^'-''^ 

^From this form one can recognize at once the Jacobi- Trudy formula [39-42] 
for the Schur polynomials S„ given as the ratio of determinants. Since Schur 
polynomials are r functions of the KP hierarchy as we had discussed in Section 
1, it is clear that 07v(A) is also r function of KP hierarchy. Due to specific 
form of this function (it contains only the odd powers of A) such tau function 
is actually tau function for the KdV hierarchy in accord with Miwa et all [34]. 
It remains to demonstrate now that Eq.s(2.17) and (2.24) also can serve as tau 
functions of KdV hierarchy. Evidently, for n ~ 2 this is the case in view of the 
arguments just presented. To prove that this is the case for n > 2 it is sufficient 
to employ the Littlewood- Richardson ("fusion" formula in physics terminology) 
rule given by 

S^.S, = ^C^,S;, (2.36) 

A 

with the Littlewood-Richardson coefficients C^^ assumed to be known [39-42] 
in principle. Successive applications of this formula to Eq.s(2.17) and (2.24) 
produces a combination of Schur polynomials each of which is tau function of 
KP hierarchy. Moreover, since for n — 2 such tau function was that for KdV 
hierarchy, evidently, for the same reasons as discussed by Di Francesco [78] the 
general case n > 2 also produces tau functions for KdV hierarchy. 
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3 Nielsen-Thurston surface automorphisms and 
partition function of 2+1 gravity 

3.1 Review of Nielsen-Thurston theory 

Let S be closed orientable Riemann surface of genus 5'. The first homotopy group, 

the fundamental group tti{S) of surface S is made of 2g generators {xi,yi}, 
i = 1 — g and a single relation so that its presentation is known to be 



Nielsen has noted that there is one to one correspondence between automor- 
phisms of 7ri(iS) and surface self-homeomorphisms. This is summarized in the 
following proposition 

Proposition 3.1. (Nielsen [68]) If g > 1, then every element of Out{TTi{S)) 
is represented by a unique isotopy class of self-homeomorphisms of S. 

An important subgroup of Out{'Ki{S)) is the mapping class group M.g dis- 
cussed in Section 1. Geometrically, this group is finitely generated by the Dehn 
twists in simple closed curves (lamination set) on S whose physical significance 
we had discussed extensively in our previous work, Ref. [4] .A simple closed curve 
C on an orientable surface S has a neighborhood £ homeomorphic to an annulus 
which is convenient to parametrize by {[7',0]|1< r < 2}. The Dehn twist in C can 
be imagined as an automorphism Tc : S S. It is given by the identity off 
£ and by [r,f?]^ [r,9 + 2nr] on £ . Using results of our previous works [1-4], 
let us illustrate these concepts on the simplest example of a punctured torus 
T2. In this case Out{TTi{T^)) = GL2(Z) and Mi,i = PSL{2,Z). Since any 
transformation from PSL{2, Z) is obtainable by projectivisation of SL{2, Z) we 
discuss everything in terms of SL{2, Z) with projectivisation at the end. Any 
transformation which belongs to SL{2, Z) is expressible in terms of 2 x 2 matrix 
A given by 



with integer coefBcients subject to condition:detA=a6 — cd=l. The character- 
istic polynomial for this matrix is given by 



This implies that the eigenvalues of A are either : 

a) both complex ( when trA=0,l,-l), 

b) both equal to ±1 (when trA=±2), 

c) distinct and real (when \trA\ >2). 

If is toral automorphism then, transformation a) is called periodic since 
(i?A)" = 1 for some n (actually, n = 12 by the Hamilton- Cay ley theorem) 



7ri(5) =< xi,yi,...,Xg,yg \ [xi,yi] • • • [xg,yg] > . 



(3.1) 




(3.2) 



- trA t -h det A = 0. 



(3.3) 
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, transformation b) is called reducible since it leaves a simple closed curve C 
invariant, transformation c) is called (pseudo) Anosov ((pseudoJAnosov if the 
line/vector field on surface (does) does not contain singularities). Such trans- 
formations are of an infinite order. Physical significance of this fact is explained 
and illustrated in our previous work, Ref.[4]. The largest of two eigenvalues is 
associated with topological entropy of the line/ vector flow and is related to the 
amount of stretching of surface <S and, hence, with the dilatation parameter of 
the Teichmiiller theory. 

Nielsen-Thurston theory generalizes the above classification of surface au- 
tomorphisms to all surfaces of genus g > 1. Already Nielsen had realized [61] 
that for 5 > 1 it is more convenient to study homeomorphisms of surface S by 
considering their image on the universal cover of <S which we choose as Poincare' 
disc model of H^, i.e. intT) U 5^ =H^.According to Nielsen [61] 

Proposition 3.2. Any lift h of the surface self- homeomorphism h: «S — > «S 
to the universal cover of S extends to a unique self-homeomorphism of the unit 
disc D, i.e. to intDUS^. 

If surface sclf-homcomorphisms h are associated with Dchn twists connected 
with a set of simple closed nonintersecting curves homotopic to geodesies (such 
set is called geodesic lamination £), then their lifts h(£) are associated with 
some maps of the circle S^^ extendable (quasi conformally) to the interior of 
the disc D as discussed in our earlier work Ref.[3] and Section 5. An image of 
the closed geodesic on S, when lifted to H^, is just a segment of a circle whose 
both ends lie on S^^ . Since geodesies are nonintersecting, circle segments on S^^ 
are also nonintersecting. In order to recover results of W-K model, in this work 
we are interested only in the periodic maps of the circle as it will be explained in 
Section 4 (after Eq.(4.4)). In connection with such maps the following remark 
is of importance. 

Remark 3.3. (A variant of Sarkovskii theorem, Ref.[79], page 88) Let /: 

^ S'"'^ be a continuous map of the circle with a periodic orbit of period 3. 
If the lift /: R ^ R has also a periodic orbit of period 3 then, / has periodic 
orbits of every period. The condition on the lift map f cannot be dropped. In 
Section 5 we argue that even though in the case of W-K model the continuous 
maps of the circle are to be replaced by the piecevise linear maps still 3 remains 
as minimal period. 

Remark 3.4. As noted by Kontsevich [9], moduli space problem makes 
sense only for Riemann surfaces obeying the following set of inequalities 

5>0,n>0,2-25-n<0 (3.4) 

with n being the number of distinct marked points (effectively distinct boundary 
components). Boundary components can be eliminated by the Schottky double 
construction. This construction can be performed as follows. If M is a com- 
plex manifold with Ci, C'n boundary components, one can consider an exact 
duplicate of it, say M, with the same number of boundary components, say. 
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Figure 3: Combinatorics of Catalan numbers presented through arrangements 
of non intersecting arcs (representing closed hyperbolic geodesies in the upper 
half plane Poincare model of H^) 



Ci,...,Cn- Evidently, for each point x G M there is a "symmetric" point x G M 
. The Schottky double 2M is formed as a disjoint union MU M and identifying 
each point x & Ci with point x & Ci for 1< i < n. In the simplest case we have 
initially either punctured torus, i.c.g = l,n = 1, or the trice punctured sphere, 
i.e. 9 = 0, n = 3. In both cases the Schottky double is a double torus. A double 
torus has 3 geodesies which belong to the geodesic lamination £. The image 
of these geodesies lifted to produces 3 circular arcs whose ends lie on 5^. 
This is minimal number of arcs required for the moduli space problem to make 
sense. According to Remark 3.3. this is also a minimal period for the periodic 
homeomorphisms of the circle in view of the Sarkovskii theorem. More on this 
topic will be discussed in Section 5. 

In the mean time we would like to discuss the general case of Riemann 
surfaces of genus g > 1 with n boundary components. It is argued in Ref.[80] 
that the total number of geodesies on the Schottky double is 65 — 6 + 3n. This 
is the dimension of space of holomorphic quadratic differentials (real on each 
of the boundary components). Hence, in accordance with Teichmiiller theory 
[20], it is the dimension of the Teichmiiller and, accordingly, the moduli space 
of such Schottky doubled surface. 

Remark 3.5. The dimension of moduli space of Schottky doubled surface 
coincides with the dimension of moduli space A^g"™** of 3-valent ribbon graphs 
used in Kontsevich paper, Ref.[9]. In the present case the same dimension for 
the moduli space (as obtained by Kontsevich) is obtained without explicit use of 
the quadratic differentials. Moreover, the equivalent of the Kontsevich-Penner 
ribbon graphs are to be obtained below in Section 4. 

It is convenient to map the circle at infinity Sl^ into the real axis R. Ac- 
cordingly, the arcs corresponding to closed geodesies on the Schottky doubled 
Riemann surface will become semicircles whose ends are located on the real axis 
R. This is depicted in Fig. 3. 

The mathematical problem associated with arrangement of the arcs dc;picted 
in Fig. 3 can be formulated according to Stanley [42] as follows: It is required to 
find a number of ways to connect 2n points in the plane lying on a horizontal 
line by n nonintersccting arcs, each arc connecting two of the points and lying 
above the points. The solution of this problem is just the Catalan number, C„ = 
/ 2n \ 

^^J^ ( j . It was used before in connection with construction of the partition 
function for our 2+1 gravity model [1,2]. Catalan numbers are very helpful in 
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Figure 4: Construction of a typical meander 



solving the mathematical problem of enumeration of meanders. Meanders had 
been used as well for description of the partition function of 2+1 gravity [1] and 
other useful statistical mechanical, dynamical and biological models [81,82]. 
Meanders can be easily constructed from the double set of arcs as illustrated in 
Fig.4 taken from our earlier work, Ref. [1] . 

3.2 Pctrtition function of 2+1 gravity 

W-K treatment of 2 dimensional topological quantum gravity is done in the 
grand canonical formalism. This means that such treatment requires two steps. 
First, one should construct the (canonical) partition function for the fixed genus 
Ricmann surface. Second, one should perform summation over all genera with 
some chemical potential. Evidently, the second step can be performed only 
if the results of the first step are available. For this reason in our previous 
works [1,2,4] only the first step was considered. Since the main Kontsevich 
identity, Eq.(1.35), is written for the canonical (fixed genus) case, it is sufficient 
to consider the fixed genus case in this work also. In order to reach an agreement 
with W-K results, it is necessary to reconsider our earlier obtained results for 
the canonical partition function of 2+1 gravity. 

Let us begin with some reminders. Meander of order n is a closed nonselfin- 
tersecting curve which intersects some straight line in exactly 2n preassigned 
points. In Ref.[l] we had discussed the way meander can be constructed from 
2 arc systems, e.g. like those depicted in Fig.3. For reader's convenience we 
reproduce a fragment of such construction in Fig.4. It is clear from this figure 
that, in general, such procedure of constructing meanders will yield a set of dis- 

(k) 

connected meanders. For a given fixed number n let the number M„ ' denote 
the total number of disconnected topologically distinct meanders whose total 
number is k. It is clear, that 1< k < n, and that 

C„ < M„ < (3.5) 

where A'/„ = Mn^~^K Each meander configuration has some statistical weight 
a;=exp{— /3J} (with /? being some fictitious inverse temperature and J is related 
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to the surface energy) dictated by the physics of the problem [1,4] so that the 
total canonical partition function Z{x) is given by 

oo n 

^5W = E^"E^'^»'^5' (3.6) 

n=0 fe=l 

with g being determined implicitly through the equation 

<k>=g-^lnZg{x) (3.7) 

with < k > denoting the average number of meanders in the cluster fc.This 
number is expected to be assigned. If this is not the case, the partition function 
should be written differently. 

Remark 3.6. The partition fimction Zg{x) is written for the system of 
meanders forming a measured foliation on a Riemann surface S of fixed genus 
g. Being guided by the Proposition 3.2., we would like to consider the lift of 
such foliation to the universal cover of S , i.e. to the unit disc model of H^. 

Remark 3.7. In both Ref.[l] and Ref.[4] the main interest in obtaining the 
partition function was to study dynamical transition from the pseudo-Anosov 
( hyperbolic) to periodic (Seifert- fibered) regime. In the present case only the 
periodic(i.e.Seifert fibered) regime is studied. This is explained in Section 4. 
Seifert fibered regime is not discussed in Refs.[l,4]. Only in this ( periodic) 
regime W-K results can he recovered from, 2+1 gravity. 

To construct the partition function on the unit disc, it is convenient to map 
the unit disc into the upper half plane model for H^. Then, the system of arcs 
representing geodesies on S is mapped into that depicted in Fig. 3. Next, the 
arc system depicted in Fig. 3 is mapped into the associated random walk as 
depicted in Fig. 5. Because of such mapping, it becomes possible to connect the 
combinatorics of vicious walkers discussed in Section 2 with that of arcs and 
meanders. Indeed, the walk depicted in Fig. 5 is directed and is not allowed to 
intersect x axis, except at initial and final points. By analogy with Fig. 4 one 
can think about another Dyck walk (below x axis). Since both walks intersect 
each other only at initial and final points this situation looks almost the same 
as for two vicious walkers. Following Labelle [83], we can make it identical to 
that for the vicious walkers. To achieve this goal, it is sufficient to translate 
the upper part by the vector (1,1) while the lower part by the vector (1,-1) so 
that one obtains either a problem about statistics of one vicious walker in the 
presence of the absorbing wall or about statistics of two vicious walkers. Both 
problems are discussed in Fisher's paper [52] and are actually equivalent. Hence, 
results of Section 2 can be applied now and one obtains the Tracy- Widom kernel, 
Eq.(2.18), in the asymptotic limit of large genus or large number of boundary 
components. This is obviously not sufficient. To go beyond this simple minded 
result requires to make several nontrivial mappings (bijections). We shall be 
brief in describing these bijections since details of the proofs can be found in 
the published literature. 
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Figure 5: A typical Dyck path [42] from (0,0) to (2n,0) made of steps (1,1) and 
(1,-1) never falling below x-axis is in one to one correspondence with a typical 
arrangement of arcs. 



We begin with the observation that to each Dyck path, e.g. that in Fig. 5, one 
can associate the Dyck word so that the Dyck path having length 2n is encoded 
by a Dyck word of length 2n. The word is composed of letters x and x in such a 
way that each North-East (respectively South-East) step correspond to the letter 
X (respectively x). The peaks (respectively troughs) correspond to the factors 
XX (respectively 5cx). Instead of North-East (respectively South-East) steps it is 
possible to choose strictly North (respectively East) steps for the Dyck paths to 
get a configuration like that depicted in Fig. 2. Hence, we obtain the following 
set of bijections: a) from the set of arcs to the set of Dyck paths; b) from set of 
Dyck paths to the set of Dyck words; c) from the set of Dyck words to the set of 
lattice paths from (0, 0) to (n, n) with steps (0, 1) or (1, 0) never rising above the 
line X = y { incidentally one can construct instead the lattice paths which never 
go below the diagonal x = y [84] ).To this set of bijections we need to add two 
more now. The first one is between the Dyck path, e.g. that depicted in Fig. 5, 
and the parallelogram polyomino. Polyomino can be made out of squares (called 
cells). A finite connected union of cells such that the interior is also connected 
and there are no cut points is called parallelogram polyomino. It is defined with 
accuracy up to translation. As depicted in Fig. 6, a parallelogram polyomino is 
bordered by two non- intersecting paths having only North and East steps. Fig. 6 
depicts bijections between the Dyck path, the Dyck word w and parallelogram 
polyomino f3{w). In such bijection the magnitude and the order of peaks and 
troughs determine the shape of the polyomino [85]. Parallelogram polyomino 
is in one to one correspondence with skew Young (or Ferres) diagram of shape 
X/ ^ . In the example displayed in Fig. 6 we have the shape (4,4,4,2,2,l)/(3,2). 
That is from standard looking Young table of shape (4,4,4,2,2,1) a piece in the 
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Figure 6: Bijections between the Dyck path, the Dyck word and parallelogram 
polyomino 



upper left corner is truncated which is also standard table of shape (3,2). Hence, 
the skew Young diagrams differ very inessentially from standard looking Young 
diagrams. Using this circumstance, we need to exhibit yet another bijection. It 
is the most important for our development. To this purpose we need to use the 
Theorem 1.11. and Corollary 1.12. in order to state yet another 

Theorem 3.8. There is a weight -preserving bijection ip between non- 
intersecting paths (Pi, ...,Pa:) and column strict Young tableaux of shape A with 
entries from N. 

To demonstrate that this is indeed the case, we follow the example discussed 
in Ref.[51] .More general case of skew column strict Young diagram is discussed 
in Ref. [42] (section 7.16). Following Ref.[51] we take fc = 4, A = (5,3,2,2) 
and N=6. Although we have now 4 nonintersecting paths their trajectories are 
completely specified by labeling of their horizontal edges (y-coordinates), e.g. 
see Fig. 2. The information about the path V^+i-i can be placed (encoded) 
into the row i of tableau T. In our case the (vicious) paths are depicted in Fig. 7 
while the corresponding encoding of the Young table is depicted in Fig. 8. Thus, 
using results of Section 1, especially Theorem 1.11 and Corollary 1.12, we obtain 
the determinant det M which contains all the configurational information about 
collection of vicious walkers. Using results of Section 2 this information is being 
translated into that for n-point correlation functions i?„, Eq.(2.24), and r„, 
Eq.(2.17). Using this information and by analogy with Eq.(3.6) the partition 
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Figure 7: Four vicious walkers encoded by their horizontal edges 
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Figure 8: The Young tableau which encodes the information about traject 
of vicious walkers 
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function can be written in the following tentative form: 



oo „ n oo 

where, the boldface for x and y variables reflects the fact that they are multi- 
dimensional, e.g. x-y=^"^j^ .Tiyj, etc. Wc use the word "tentative" because we 
would like at this point to make a connection with recent works by Okounkov 
[77,86] and Okounkov and Pandharipande [87] who use the asymptotics of Hur- 
witz numbers to arrive at results similar to our Eq. (3.8). For reasons which 
will become clear in the next section and for the sake of agreement with these 
recent works, Eq.(3.8) should be rewritten as follows: 

Z„[x] =Z[xi,...,x„] = ^-4^^£^ ) (3.9) 

i 

with 

/n 
[]d|/,exp(x-y)T„(y). (3.10) 
i 

Remark 3.9. For developments presented in this paper Eq.s (3.9) and (3.10) 
are actually unnecessary. They are provided, nevertheless, because they are 
interesting in their own right and, although the initial arguments in these papers 
are different from ours, the end results, e.g. Eq.(3.10), is in accord with our 
Eq.(3.8). 

The interest in these equations stems from the fact that in the large n limit, 
Eq.(3.10) acquires familiar in physics literature path integral form. This fact 
may lead to some new physical applications. Already many physical applications 
of vicious walker models can be found in papers by Fisher [52] and by Huse 
and Fisher [53]. 

In the present case, using results of Tracy and Widom ( Eq.(4.6) of Ref.[56]) 
Eq.(2.18) of Section 2 can be rewritten as follows 

oo 

K{X, y) = j dzAi{X + z)Ai{z + Y) (3.11) 



with Airy function defined in Eq.(2.19). Using Okounkov's Lemma 2.6, Ref. 
[77], 



OO 

/ 



, X-/ ^^•/ ,x / N 1 a + b (a-&)^, 
dzAi{z + a)Ai{z + 6) e^^{xz) = —= exp(— —x ), 

(3.12) 
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Eq.(3.10) can be rewritten as follows: 

g(x)= K^^^"^ /nc^y.exp(-Vi^i-^-V^i±^.0 
^ ' 2"7rt Yl^i y 1.1 Z_. 4a;. ^2 

(3.13) 

provided that the cyclic boundary condition , = j/i, is imposed. This makes 
the above "path integral" reminiscent to that for the ring polymer near hard 
nonpenetrable wall in the presence of some stretching force [67]. 

To make connections with results of Kontsevich [9] we need to rewrite the 
partition function Z[x] as genus expansion: 

oo 

Z[x]=^Z,(x) (3.14) 

5=0 

where, in view of Eq.(1.2), 

n 

Zg{^)= E <rd,---rd^>gj{xf. (3.15) 

E(ii=3g-3+n i=l 

The question arises immediately about connections of Eq.s (3.14) and (3.15) 
with Eq.s (3.8)-(3.10). This issue is addressed and treated in Refs. [77,87] and, 
therefore, there is no need to repeat the arguments presented in these references 
here. Instead, to reach essentially the same goals we would like to use different 
arguments in this work. 

In Section 1.1. wc had noticed that the calculation of tautological averages is 
effectively equivalent to the calculation of the Weil-Petersson volumes of ■Mg,n- 
In algebraic geometry there is a Wirtinger formula for such volume calculations 
[12]. In connection with Grassmannians and Schubert calculus this formula 
had been discussed in the fundamental paper by Chern [46]. Wirtinger- like 
formula is used in Kontsevich paper as well (e.g. see Section 3 of Ref.[9]).Hence, 
the partition function of two-dimensional topological gravity is effectively the 
generating function for the Weil-Pctersson volumes. In the fixed genus case such 
volume according to Kontsevich is given by 

volp^{Mg,n) = J {plci{Xi) + ---+plc{Xn)f (3.16a) 

with ci(Aj) being the first Chern class of the i— th line bundle, i = 1, ...,n and 
d = 3g — + n. The above expression becomes a true volume when in the set p* 
of indeterminates p* ={pi, ...,Pn) (arbitrary sequence of positive numbers) each 
Pi is being put equal to one. This is not necessary, however, since the quantity 
of interest is the product given by Eq.(1.2) which is obtainable anyway with 
PjS being different from one. Hence, the indeterminates actually play a role of 
an auxiliary variables analogous to Xi in Eq.(3.15). The connection with Schur 



30 



polynomials and, hence, with tau function is clear if one combines Eq.(1.28) and 
Porteous formula, Eq.(1.27) with Eq.(3.16a). Finally, arguments presented in 
Section 2, especially, Eq.(2.36), provide needed justification of Eq.(3.15) since 
the result, Eq.(3.16), of Kontsevich can be equivalently rewritten as 

volp*{Mg,„) = sgnx ^ < ^d.i-Td^ > ■ (3.16b) 

E<ii=<i i=l 

To make connection with the matrix models additional steps are required. For 
instance, Kontsevich is making a Laplace transform of Eq.(3.16b) in order to 

obtain 

C{volp,{Mg,nmi,...,Xn) = E < ^<i^-^<ir^ > n^Ar^^'^' + l) ^3_^7^ 

Sdi=d i=l 

with Aj being the Laplace variable conjugate to pi. Taking into account that 
{2di)\/di\ = 2"^* ■ {2di — 1)!! the overall factor of 2** drops out from the main 
Kontsevich identity, Eq.(1.35). In view of the results of Refs. [77,86,87] and that 
going to be presented in Section 4, in order to reach an accord with the result 
of Kontsevich, Eq.(3.17), it is also necessary to perform the Laplace transform 
on ^g(x) in Eq.(3.15) provided that this expression is properly rescaled. The 
Laplace transform is obtained then as follows: 

°° n 

£(Zg(x))(a,...,U= /TTrfa;iexp(-^x)Zg(x) -. (3.18) 

ill (27r)^nr=i^f 

Application of the formula 

oo 

J dyy--^ eM-sy) = r(n + i)."""^ = V^ ^f^"^'! (3.19) 



to Eq.(3.18) produces the expected result: 

£(Z,(x))(a,...,U= E <^ci.-.-r,„>n ^^^'~J-\ (3.20) 

where, in order to achieve an agreement with Kontsevich, one needs to make 
identification: Ai = \/2^ in Eq.(1.35). Clearly, such identification is ultimately 
connected with rescaling made in Eq. (3. 18). The justification of this rcscaling is 
explained in the next section. This is also needed for completion of the proof of 
Kontsevich identity, Eq.(1.35) of Section 1. 



4 Ribbon graphs,Young tableaux and Kontse- 
vich identity 
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4.1 Construction of ribbon graphs 

Ribbon graphs had been invented by Penner [ 62] ( similar construction can 
be found also in papers by Harer [88] ) for description of the moduli space 
■Mg^n of Riemann surfaces. In physics literature similar construction had been 
independently developed by Saadi and Zwicbach [63] and later developed in nu- 
merous papers by Zwiebach, e.g. see Ref. [89] cited in Kontsevich work, Ref.[9]. 
Both Zwiebach and Kontsevich use the Jenkins-Strebel quadratic differentials 
(e.g. sec Section 3 of our earlier work, Ref. [2], for condensed summary of their 
properties) for construction of the ribbon graphs. In this section we would like 
to develop alternative method of construction of ribbon graphs which does not 
require explicit use of quadratic differentials. We would like to explain the ratio- 
nale for constructing the ribbon graphs in connection with results obtained thus 
far in this paper. First, the ribbon graphs appear naturally in matrix models 
for strings [90] and QCD [91]. Hence, they arc physically interesting and easily 
constructible using Feynman-like rules known in quantum field theory [92] . Sec- 
ond, mathematically these graphs are interesting for several reasons : a) they 
arc associated with problem of imbedding of, say, trivalent graphs into Riemann 
surface of fixed genus [93] and, since graphs are related to combinatorial group 
theory [94], this problem is also of group-theoretic interest; b) these graphs 
are also related to train tracks discussed in our work, Ref. [2]. Because of this 
connection with train tracks, the number of potential practical applications is 
expected to be well beyond particle physics as explained in the discussion sec- 
tion of Ref. [2]. Third, and most important for the purposes of this work, they 
are needed for establishing one-to one correspondence between the topological 
gravity and the random matrix models of string theory. The aspects of this 
correspondence are discussed in detail in Ref. [90]. Mathematically, this corre- 
spondence is reflected in the Kontsevich identity, Eq.(1.35). Hence, following 
Kontsevich [9] it is necessary to prove that the r.h.s.of Eq.(1.35) is equal to the 
l.h.s. In sections 2 and 3 we had demonstrated that the l.h.s of the identity, 
Eq.(1.35), is associated with the enumeration of allowed configurations for an 
assembly of vicious walkers. This problem has been mapped into enumeration 
of the Yoimg tableau, e.g. see Figs 7 and 8. These pictures provide a geomet- 
rical way of describing partitions of nonnegative integers. Hence, for ribbon 
graphs we have to find analogous partitions. Evidently, the Kontsevich identity, 
Eq.(1.35), is just the statement about existence of different but equivalent ways 
of describing partitions as had been stated in the Introduction. 

To describe partitions associated with ribbon graphs we need to construct 
such graphs first. To facilitate reader's understanding, we employ some results 
on graphs and ribbon graphs from the pedagogically written paper by Mulase 
and Penkava [64] .Ultimately,our way of constructing the ribbon graphs is dif- 
ferent from that discussed in this reference since we are not using quadratic 
differentials. 

Definition 4.1. A graph r=(V, £ ; i) consists of a finite set of vertices 
V=( Vi,..., Vm) and finite set of edges S together with a map i from £ to the set 
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(V X V) /S2 of unordered pairs of vertices called incidence relation. The quantity 

aij = \i-\Vi,Vj)\ (4.1) 

is the number of edges connecting vertices Vi and Vj . The degree (valence) of 
a vertex Vj is the number 

deg{Vj)=Y,ajk + 2ajj (4.2) 

which is the number of edges incident to the vertex. 

Definition 4.2. A graph isomorphism is a pair (a,/3) of bijective maps 
a : V ^ V'and (} : £ ^ £' that preserve the incidence relation. 

To construct a ribbon graph the above definitions should be modified. The 
modification consists in labeling of middle points of each edge thus effectively 
creating extra degree 2 vertices associated with each edge. We denote this extra 
vertex set as Vg. Now the new set of vertices is the disjoint union V U while 
the new set of edges is the disjoint union £ U£ since the midpoint of each edge 
now divides it into two parts. The incidence relation is described now by the 
map 

is:£u£ — ^ V X Vf (4.3) 

because each edge of is connecting one vertex of V to one vertex of Vs 
.Obviously, an edge of Ts is called a half edge of F. For every vertex V& V of F 
, the set i£^({l^} x Vs) consists of half edges incident to V so that 

deg{V) = %\{V}xVs)\ (4.4) 

Definition 4.3. A ribbon graph is a graph F=(V, £] i) together with a 

cyclic ordering on the set of half-edges as depicted in Fig. 9. 

Such definition of the ribbon graph leads to the following construction of 
such graph out of ordered vertices. The strips corresponding to the two half 
edges arc connected following the orientation of their boundaries to form rib- 
bons. The final surface is no longer planar in general. It is an oriented surface 
whose boundaries are made of boundaries of the ribbons. To illustrate this 
construction, let us consider the simplest case first. To this purpose, we re- 
call that for the case of punctured torus and trice punctured sphere we have 3 
geodesies (on the Schottky doubled surface) which we lift to the Poincare disc. 
We associate with these geodesies a trivalent graph as depicted in Fig. 10. 

Next, we make another copy of this picture. Then, we thicken the edges 
emanating from the vertices and provide orientation in accord with Fig. 9. Next, 
we glue the strips to each other. This can be done in two ways. One is depicted 
in Fig. 11 while the other is depicted in Fig. 12. We would like to notice that 
in both cases we have constructed ribbon graphs Fg,„ = F with the number 
of thickened edges equal to 6.g — 6 + 3n number of edges (e.g. 3 for both the 
trice punctured sphere and the punctured torus) in accord with Kontsevich 
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Figure 9: A cyclic ordering at the vertex of the ribbon graph as compared with 
no ordering in the case of a vertex for ordinary graph 



Figure 10: An elementary 3 valent vertex associated with geodesies lifted to the 
Poincare disc 
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Figure 11: One way to make the simplest ribbion graph (the underlying "nor- 
mal" graph is shown for comparison) which is topologically equivalent to the 
punctured torus. Notice, that the same underlying "normal" graph also gives 
rise to another ribbon graph which is topologically equivalent to the trice punc- 
tured sphere (not to be confused with that depicted in Fig. 12). 



Figure 12: Another way of making the simplest ribbon graph (the underlying 
"normal" graph is also shown for comparison). The ribbon graph is topologically 
eqivalent to the trice punctured sphere 
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Figure 13: Typical contraction-expansion (Whitehead) moves characteristic for 
W-K model 



[9]. To obtain more complicated graphs it is convenient to proceed by induction. 
To this purpose, following again Ref.[64], we need to define the operations of 
contraction and expansion. 

Definition 4.4. If the edge E of F is incident to two distinct vertices Vi 
and V2 another ribbon graph T' called contraction of F is obtained from F by 
removing the edge E and joining the vertices Vi and V2 to a single vertex with 
the cyclic ordering at the joint vertex determined by the cyclic order of the 
edges incident to Vi starting with the edge following E up to the edge preceding 
E, followed by the edges incident to V2 starting with the edge following E and 
ending with the edge preceding E. 

Evidently, the contraction procedure decreases the number of edges and 
vertices by one. Every ribbon graph can be obtained from the trivalent ribbon 
graph by applying a sequence of contractions. 

Definition 4.5. The expansion of the ribbon graph is operation inverse to 
the contraction. 

This means that every trivalent ribbon graph can be constructed by the 
expansion procedure starting from a very simple analogue of Fig. 10. In view of 
the Euler equation 

V(F)-£:(F) = 2-2g-n (4.5) 

such procedure will produce the desired Kontsevich-type trivalent ribbon graph. 
Let us consider this construction in some detail. For example some representa- 
tive contraction and expansion of an edge is depicted in Fig. 13 . 

Remark 4.6. The reader familiar with our earlier work, Ref.[l], can easily 
recognize the Whitehead moves characteristic for train tracks. This fact alone 
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Figure 14: Expansion/contraction of a representative vertex facilitated by the 
associated with it dual polygon 



is sufficient for making connection between ribbon graphs and quadratic differ- 
entials. Since in the present case no topology changing moves are involved, this 
justifies our earlier statement that the results of W-K can be recovered from the 
Seifert-fibered (periodic) phase of 2+1 gravity. 

For any vertex of degree > 4 there are d{d — 3) /2 ways of expanding it 
by adding an edge. Consider a portion of T made of vertex of degree d > 4 
and d half edges emanating from this vertex. A dual to this portion is a convex 
polygon of d sides depicted in Fig. 14 along with the underlying vertex. The 
series of contractions/expansions is ultimately related to the number of ways a 
convex polygon with d sides can be triangulated by nonintersecting diagonals. 
This number is Catalan number again [42]. Although one may probably use 
this observation to develop bijections analogous to those discussed in Section 3, 
we leave such an opportunity outside the scope of this work in view of another 
options to be discussed below. 

To finish our construction of the ribbon graphs we have to keep in mind 
few additional facts. First, by looking at Fig. 10 and realizing that two copies 
of the disc model of are needed for construction of the ribbon graph, it is 
clear that the process of expansion and gluing of strips to each other should be 
allowed to be different in each disc. Hence, we have to associate the probability 
of 1/2 to every vertex on each disc. This means that we have to have an overall 
(assembly) factor of 2^^'-'") for a particular realization of the ribbon graph. 
Second, the rest of the automorphisms of the ribbon graph are, evidently, the 
same as for the "normal" graph. 

This concludes our description of the ribbon graph construction. 

4.2 Young tableaux and Kontsevich identity 

In section 3 the Laplace transform, Eq.(3.20), of the partition function for an 
assembly of vicious walkers has been obtained. Now it is time to explain how 
these vicious walkers reemerge with help of the ribbon graphs. To this purpose 
let us notice that the boundary components of these graphs can be looked upon 
as made of polygons as Kontsevich had noted in section 2.2. of Ref.[9]. Hence, 
each ribbon graph is can be classified by certain non negative integer number of 
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triangles (na), quadrangles (n4), pentagons (ns), hexagons (ng), etc. Each edge 
Ei{T) of the ribbon graph has some length li. Evidently, these lengths can be 
grouped into sets identified with polygons so that the numbers above represent 
the multiplicities for these sets. Suppose, that the total length L of all polygons 
is prescribed in advance. Then, given ribbon graph can be viewed as particular 
realization of the partition of L into numbers associated with lengths of these 
polygons. Alternatively, one can assign the total number of polygons (that is 
the total number of faces n in Eq.(4.5)) and consider partition of this number 
into ns, ns, ne , etc. In any case, this means that the Young tableaux can be 
associated with such partition and, in view of the results of Section 3, again, 
the vicious walkers can be linked with such Young/Ferrers tables. The question 
remains: how to connect the vicious walkers with ribbon polygons explicitly? 
To this purpose recall Remark 2.1. and discussion which follows. According 
to this remark and the following discussion both correlation functions i?„ and 
T„ are defined in fact for vicious walks made out of loops. Topologically, our 
polygons are also loops. Hence, it makes sense to identify the loops entering 
into expressions for i?„ and T„ with those coming from polygons associated 
with ribbon graphs. This identification requires some care (that is it requires 
some proofs) and, hence, cannot be made straightforwardly as we would like to 
demonstrate now. 

If we associate (replace) the polygonal paths by the paths of vicious walkers 

then, except at vertices, only "binary interactions" between these walkers need 
to be considered. These "binary interactions" are of geometrical origin since 
they are effectively equivalent to considering just one vicious walker in the pres- 
ence of an absorbing wall [52]. According to Fisher [52], this means that "no 
walk can penetrate the wall and any walk attempting to do so is eliminated" . 
In addition to the geometrical constraint on such walk, the present case differs 
from that discussed by Fisher because one has to demonstrate that the walk 
which survived encounter with one wall entirely forgets about this encounter 
when it is facing the next wall. Fortunately, this is the case. Indeed, the non 
normalized distribution function Qn{x:a) for the walk starting at x = a and 
after n steps reaching the point x is obtained by Fisher, e.g. see Eq.(5.5) of 
Ref. [52] , and is given by 



where, apart from the factor aexp{—an),oiie can easily recognize the probabil- 
ity of the first passage through x at "time" n for the random walk which had 
started at the origin: x = [75]. This probability has yet another interpreta- 
tion more suitable for the problem we are discussing. Indeed, following Feller, 
Ref.[75],Chr3, the same expression describes probability for two dimensional 
(directed) random walkers which had begun their walk at the origin and af- 
ter n steps had ended at some point a: > 0, y > of a;, y- plane, provided that 
these walkers never cross the x— axis. Evidently, one may associate a;— axis with 
"time" (in our case n) direction while y -axis with "space" direction in order 



Qn{x;(T) 



e-'"" aa:exp(-f^) 



(4.6) 
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to get vicious walk interpretation of such random walk. In any event, it should 
be clear that in space direction (that is perpendicular to the wall) such walk is 
acting as random and, hence, it "forgets" its past. This observation provides 
justification for representation of the "circular" vicious walk by a sequence of 
independent random vicious walkers- one for each edge of the closed polygon. 

Next, we have to find an analogue of the partition function, Eq. (3. 8). This is 
accomplished in several steps. First, wc notice that Eq.(3.8) docs not contain 
information about the number of steps in the walk, only about the number of 
walks. Using Eq.s (3.8) and (4.6) we introduce for each edge the following 
generating function (known as perimeter measure [86,87] ) 

oo 

C(a,6)= J dxP{^,,^2;x) (4.7) 



where 

oo oo 

-P(Ci,6;a;) = y dm j dn2Qm{x;^i)Qn2{x;C2) (4.8) 



and we put a = 1 in Qn{x;0 Eq.(4.6). This is permissible since the factor 

exp{—an) can be always amended to absorb a (e.g. read Fisher's paper, Ref.[ 
] ). Using Eqs (4.6)-(4.8) the explicit form of the perimeter measure can be 
easily obtained using standard tables of the Laplace transform. The final result 
is given by 

C(a,6)= (4.9) 

This result should be applied to all edges of the ribbon graph and the final result 
should take into account all automorphisms |Azttr| of the underlying "normal" 
trivalent graphs. Taking into account the assembly factor of 2"^^'"' the final 
result can be written as follows 

2-v{r)2S(r) i 
£(Z,(x))(a,...,^„)= E \Aut(G)\ n /Of , /o, (4-10) 

with indices i and j referring to ■i-th and j— th polygons sharing the same edge. 

Since the l.h.s. is given by Eq.(3.20), Eq.(4.10) coincides with Kontsevich 
identity, Eq.(1.35), provided that identification = is made. Since for 

trivalent graphs |V| = | |f | , we obtain 

2-v(r)2£(r) ^ 2'^3-2+n ^^^^^^ 

in accord with Ref.[87]. 
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5 Piecewise linear homeomorphisms of the cir- 
cle and KdV 



In this section we would like to provide another interpretation to the results 
obtained in previous sections. This interpretation is desirable since its aim is 
to explain to what extent W-K model is universal from the point of view of 
dynamical systems theory. This is also desirable for the purpose of connecting 
results of this paper with WDVV equations and Frobenius manifolds [27,95]. 
To begin, we need to introduce some information about the Thompson groups. 

5.1 Some facts about the Thompson groups 

The groups F, T and V were introduced by Richard Thompson in 1965. Unfor- 
tunately, his results had not been published. This, nevertheless, had not stopped 
line of research initiated by Thompson as can be seen from the pedagogically 
written review article by Cannon, Floyd and Parry [96]. Below, we provide 
some basic facts on Thompson groups using mainly this reference. Additional 
sources will be used whenever they are needed to suit our needs. 

Let F be the set of piecewise linear homeomorphisms from the closed unit 
interval [0,1] to itself that are differentiable except at finitely many points. Let 
f G F and let = a;o < a^i • •• < a;„ = 1 be the set of points at which / is not 
differentiable. This partition determines intervals [xi_i, Xi] for i = 1, n which 
are called intervals of the partition. A partition of [0,1] is called standard dyadic 
partition if and only if the intervals of partition are standard dyadic intervals. 

Definition 5.1. A standard dyadic interval in [0,1] is an interval of the form 
[^5 where a and n are nonnegative integers with a < 2" — 1. 

It is useful to associate a finite ordered rooted binary tree associated with 
standard dyadic intervals as depicted in Fig. 15. For n, 
the function / can be written as follows 

f{x)=aiX + bi (5.1) 

with Ui being a power of 2 and bi being a dyadic rational number. It can be 
shown that e F and / maps the set of dyadic rational numbers bijectively 
to itself. Hence, F is closed under the composition of functions and therefore is 
a subgroup of the group of all homeomorphisms from [0,1] to [0,1]. This group 
is F group of Thompson. 

Definition 5.2. When points and 1 are identified to make a circle S^, 
then, the resulting Thompson group is called T group. Another Thompson 
group acting on the circle is V group. Its definition is a bit technical [96] but, 
at "physical level" of rigor, the difference between V and T groups is hardly 
noticeable. Therefore, following Ref.[97] we shall denote both V and T groups as 
PL2(S'^) and we shall keep in mind that both groups are subgroups of the group 
Homeo+ of piecewise orientation-preserving homeomorphisms of the circle. The 
following proposition proven in Ref.[96] is very important. 
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Figure 15: The tree of standard dyadic intervals 



Proposition 5.3. Let / G F. Then there exists a standard dyadic partition 
= < xi • •• < a;„ = 1 such that / is Unear on every interval of the partition 
and = f{xo) < f{x\) ■ ■■ < f{xn) = 1 is a standard dyadic partition. 

Evidently, using the proof of this proposition, the same statements can be 
made for Thompson groups T and V . 

Remark 5. 4. Using results of Ref.[42], the tree depicted in Fig. 15 can be 
put in bijective correspondence with the set of non intersecting arcs depicted 
in Fig. 3. In view of Proposition 5.3., different arc configurations correspond 
to different piecewise linear homeomorphisms of the circle caused by PL2(S'^). 
This statement is going to be examined further below in Section 5.3. 

Remark 5.5. In view of Eq.(5.1) it can be shown [98] that, except for points 
= a;o < xi • •• < a;„ = 1, the group PL2(S'^) is isomorphic to group PSL2(Z) 
which had been discussed at length in our previous work, Ref.[4]. 

5.2 Some facts about the Ptolemy group 

In addition to isomorphism mentioned in Remark 5.5. there is yet another 
important isomorphism between the Thompson and the Ptolemy groups which 
we would like to describe briefly in this subsection. Following Lochak and 
Schneps [99] define the standard marked tessellation of the Poincare disc as 
dyadic tessellation with marked (oriented) edge from to oo as depicted in 
Fig.16. 

The elementary move a on the oriented edge of the dyadic tessellation 
changes its location from one diagonal of the unique quadrilateral containing 
it to another one by turning it coimtcrclockwise as depicted in Fig. 17. Such 
elementary move is of order 4 evidently. In addition to this move there is an 



41 



Figure 16: The standard dyadic tesselation of the Poincare disc with oriented 
marl^ed edge being displayed 




Figure 17: The elementary move on the oriented edge 
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Figure 18: An arrow- moving move 



arrow-moving move /3 depicted in Fig. 18. This operation moves an oriented edge 

to another edge without changing the tessellation itself. The minimal order of 
this operation is 3. This fact makes it formally similar to that coming from 
Sarkovskii theorem, e.g. see Remark 3.3, which forbids periodic orbits of period 
lesser than 3. Generators a and (3 along with relations 

a* = l;/33 = l;(a/3f = 1 (5.2) 

and commutator relations 

[/3a/3, a^l3a[3a^] = 1; [/3a/3, Pa^ Pafia^ 13^ a^] = 1 (5.3) 

define the Ptolemy group [99].Imbert had demonstrated [97] that this group is 
isomorphic to PL2(S'^). The same result is obtained in Ref.[99] with help of 
slightly different methods. 

5.3 Circular homeomorphisms, combinatorics of the rib- 
bon graphs and volumes of hyperbolic polyhedra 

Now we would like to add some important details to Remarks 5. 4. and 5.5 made 
above. In view of the Remark 5.5. the standard dyadic tessellation of the 
Poincare disc depicted in Fig. 16 can be shown [97,98] to be in one to one corre- 
spondence with the Farey tessellation of the unit disc D discussed in detail in 
Ref.[4]. Such tessellation can be obtained from two "seeds", that is from two 
triangles with vertices located at 0=1/0, 1=1/1 , oo = 1/0 and 0,-1 and oo re- 
spectively (e.g. sec Fig. 16 and Section 2 of Ref.[4]) The duals to the tessellations 
originating from these seeds are rooted binary trees whose roots are connected 
to each other as depicted in Fig. 9 of Ref.[4]. The Farey tessellation is also known 
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as ideMi triangulation of D [98] . Ideal since the triangle which cover the hyper- 
bolic plane/disc are ideal in the same sense as Euclidean equilateral triangles 
covering the entire plane without gaps form an ideal triangulation of the Eu- 
clidean plane. Since the Farcy tessellation is multileveled (e.g. Fig. 16 depicts 
just two levels while Section 2 of Ref. [4] describes how to construct tessellations 
of any level). If the level is finite(or at least periodic as depicted in Fig. 8 of 
Ref. [4]), then the duals of these tessellations are two rooted binary trees that 
are also finite. With help of Ref. [42] they can be put in bijective correspondence 
with two systems of arcs, e.g. those depicted in Figs. 3, and 10, etc. needed for 
construction of the finite ribbon graph as discussed in Section 4. Surely, the 
infinite level case of Farey tessellation is obtained as limiting case of finite level 
case. The moves depicted in Fig.s 17 and 18 are in fact homeomorphisms of 
the circle [98,100] belonging to the subgroup Honieo^. They affect the pattern 
of triangulation of the unit disc and, hence, affect the binary trees. As demon- 
strated in Ref. [96], the group of automorphisms of the reduced binary trees ( 
meaning of the word reduced is explained in Section 2 of Ref. [96].) coincides 
with Thompson group PL2(>S'^). In the same reference the analogue of moves 
a and /3 for trees is obtained. Some additional physics can be associated now 
with these mathematical facts. To this purpose, following Ref. [101] we notice 
that the automorphisms of two binary trees associated with two arc systems 
are independent. Each of these trees is in bijective correspondence with some 
triangulation of the convex polygon by non crossing diagonals [42] . Every auto- 
morphism of individual tree causes some change in the triangulation pattern for 
given polygon. Since these changes are independent, following Sleator, Trajan 
and Thurston [102 ] one may ask about the minimal number of automorphic 
steps needed to bring one triangulation pattern in exact correspondence with 
another. This problem is of major importance in computer science. To make it 
also of importance to physics (and also to evolutionary genetics [103], dynam- 
ics and thermodynamics of folding of RNA [82,104], etc.) following Ref. [102 ] 
several more steps are required. First, by gluing these two polygons to each 
other we form a triangulated sphere S'^. Second, we project this sphere stere- 
ographically into the complex plane so that pattern of triangulations on is 
transferred to the Euclidean plane. Next, we associate this plane with boundary 
of the upper half space model for H^.The triangular pattern on such plane is 
sufficient for reconstruction of the hyperbolic polyhedra according to Thurston 
[16]. A pedagogical account of how it can be actually done can be found in 
Ref. [105 ]. Finally, if needed, such polyhedra can be mapped into the hyper- 
bolic ball model of H"^ so that the vertices of these hyperbolic polyhedra are 
located at the sphere at infinity 5^. Each triangulation pattern on the sphere 
S"^ thus associated with some hyperbolic polyhedron of finite volume Vh which 
is determined by this triangulation pattern. The physics can be now injected 
into this picture by introducing some Boltzmann factor exp(-/3Viy) so that dif- 
ferent volumes are related to different triangulations. The larger the volume 
is for a given "temperature"/? , the more stable triangulation pattern becomes. 
Incidentally, Ref. [102 ] is devoted to finding hyperbolic polyhedra with large 
volumes. 
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Remark 5.6. Very recently statistical mechanics of various physical sys- 
tems with Boltzmann factor containing volume (including hyperbolic volume) 
had been considered by Atiyah and Sutcliffe [106 ] . Actually, these authors had 
considered instead of volume its logarithm. 



5.4 Universal Teichmiiller space, KdV and PVobenius man- 
ifolds 

Connections between the Dyck paths, Virasoro characters and exactly solvable 
lattice models had been known for some time [107 ]. Moreover, since combina- 
torics of the Dyck paths is directly associated with that for the Catalan numbers, 
such connections are actually not totally unexpected since there are about 150 
interpretations of the Catalan numbers [42]. This, by the way. means that 
potentially there are much more applications of the W-K model than we had 
mentioned so far. It is not the purpose of this subsection to provide the list of 
such applications. Instead, we are interested in intrinsic features of W-K model 
which are encoded by the Thompson (or the Ptolemy) group. Following Penner 
[98,100], the universal Teichmiiller space T^ss is defined by 

Homeo+ 

Mob ^ ' 

where Mob is just PSL{2, R).li T^^g' is a set of all tessellations of D, then 
%ss — "^ss' as well, as Penner shows. In Section 7 of Ref.[3] we had 
introduced and discussed the universal Teichmiiller space T(l) as defined by 
Bers. More explicitly, such space is defined by 

with QS being a set of all quasisymmetric deformations of the circle S^^ (or 
real line R).Wc shall discuss these deformations in some detail below. In the 
meantime we notice that "physical" definition of the Teichmiiller space T is 
given by 

(5.6) 



PSL{2, R) 



where G= Diff+ denote set of all orientation preserving diffeomorphisms of 
5^. Since G is proper subgroup of QS [98] , it is clear that T is embedded into 
T(l). In Ref.[98] Penner argues that QS is subset of Homeo+ . This leads to 
the following inclusions: 

TcT(l)cTe,,. (5.7) 

This result is very nontrivial since it makes sense out of Sarkovskii theorem 
(Remark 3.3 and Section 5.2) in the present context. It also allows to use Nag 
and Verjovsky [108] arguments (summarized in Section 7 of Ref.[3]) for use 
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of QS deformations of S}^ in order to obtain the Virasoro algebra. It is well 
documented fact that the method of coadjoint orbits [109] directly connects the 
Virasoro algebra with KdV equation. The connection between QS deformations 
and KdV can be established directly without method of coadjoint orbits. For 
this purpose we need to use some classical results of Ahlfors [65] . 

Let w = f{z) be a homeomorphism of the complex ^-plane (or 5^). Then, 

dw = fzdz + fzdz. (5.8) 

The complex dilatation factor fi giving rise to the Beltrami equation is defined 
now as 

h = m/z (5.9) 

and the associated modulus of this dilatation factor is defined by rf/ = < 
l.The mapping is considered to be K-quasiconformal if there is a nonnegative 
constant K such that Df < K ,where -D/ = jz-^- Accordingly, the mapping is 
conformal if Df = 1. Suppose now that / maps the upper Poincare half plane 
to itself (which is, of course, equivalent to the mapping of D to itself). The 
mapping is quasisymmetric, i.e./ G QS if for all points x, x — t and x + t on 
real line R the following M condition 

M-i < {S^+*^~ ^^""l < M (5.10) 
- fix) -f{x-t)- ^ ' 

is satisfied with M being some non negative constant. In Chr.6 of Ref.[65 ] 
Ahlfors proves that the function which is schlicht and has a quasiconformal ex- 
tension to the upper half plane must obey the following Fuchsian-type equation 
known already to Poincare [110] 

y" + ^(Py = 0. (5.11) 

The mapping function / = yi/y2 where yi and 2/2 are being two independent 
solutions of Eq.(5.11) normalized by 

y'm - 2/22/1 = 1- (5-12) 

The function </> is determined by equation [f]=(j) with [/] being the Schwarzian 
derivative of /. Recall [111], that the function is considered to be schlicht (or 
simple) at some point of complex plane if its first derivative is nonzero at this 
point. 

Lazutkin and Pankratova [112 ] studied Eq.(5.11) from the point of view of 
the circle maps. They used a variant of Eq.(5.1) given by 

+ 2tt) = F(0 + 27r, e e R, (5.13) 

to study the transformational properties of Eq.(5.11). In particular, change of 
variables 



F(e) andy(i^(0)=r(0/nO (5.14) 
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leaves Eq.(5.11) in the same form Y" + = with potential $ given by 

m = <)>iF{0)[F'iof+[m] (5.15) 

with again, being the Scliwarzian derivative of F. Eq.(5.15) actually de- 

termines transformational properties of the Seliwarzian derivative. In Ahlfors 
notations [65] this equation can be equivalently rewritten as equation for trans- 
formation of the Schwarzian derivative: 

[foF] = {[f]oF]){F'f + [F]. (5.16) 

Evidently, this equation holds irrespective to the explicit form of Eq. (5. 13). Using 
this observation it is of interest to consider transformation of the type i^(^) = 
^ + 6(fi{^) with S being small parameter and <^(^) some function which makes 
Eq.(5.13) to hold. Use of this type of function in Eq.(5.16) produces 



*(C + S^iO) = m + S [T^) (0 + 0{d') (5.17) 

with 

(r.^) (0 = m^'iO + + cf^'iO^iO- (5.18) 

The Korteweg- de Vries equation can be obtained now in a simple minded way. 

To this purpose, since ip is arbitrary, we choose = <f>. Next, we make both ip 
and (/> to depend upon the parameter 5 , that is we write (j) = (/>(^, (5). Next, we 
assume that the parameter 5 plays role of "time" t and, finally, we write 

lim ^ = 5.19 

t^o t dt ^ ' 

This then produces our final result: 

^=3<l>cP'+^-r- (5.20) 

Eq.(5.20) coincides with Eq.(2.2) of Segal [113] and, hence, can be called KdV 
equation. We provided details of derivation in order to emphasize the universal- 
ity of this equation in problems which involve circular maps or maps of D. Since 
KdV is effectively dual to the Virasoro algebra surely it also can be obtained 
via Nag-Verjovsky approach to construction of the Virasoro algebra and the 
Kirillov-Kostant two- form by using the universal Teichmiiller space T(l) [108 
] . Summary of Nag-Verjovsky results can be found in our earlier work, Ref. [3]. 
These results were extended by Penner [98] whose Theorem 5.5. essentially as- 
sures the existence of the Kirillov-Kostant two form which is invariant under 
transformation taken from the Ptolemy group. Independent and deep studies of 
the same problem can be found also in earlier work by Ghys and Sergiescu [114 
] . Higher order KdV can be easily obtained by sequential use of T operator as 
can be seen from Ref. [115 ]. 
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Eq.(5.11) contains actually much more information than wc had discussed 
thus far as was shown by Dubrovin, e.g. see Example C.l.of Ref.[95]. From 
this example (and preceeding discussion) it follows that variety of equations 
of Wittcn-Dijkgraaf-Vcrlindc-Vcrlindc (WDVV) type arc obtainable as special 
cases of Eq.(5.11). Moreover, the self-dual Yang-Mills and Einstein equations 
also follow from Eq.(5.11) [95,116]. Geometrically, solutions of WDW equa- 
tions represent Frobenius maifolds. Thy arc named after Frobcnius who was 
the first to discover them in 1882 [117]. Incidentally, the self-dual Einstein 
equations had been studied already in 1881 by Halphen [118] and rediscovered 
by Atiyah and Hitchin [116]. In recent paper by van dc Lcur and Martini [119] 
KP representation theory and the related Sato infinite Grassmanian are used 
to construct solutions of the WDW equations and, hence, the Frobenius mani- 
folds. In addition, Dubrovin and Zhang [120] had recently demonstrated that , 
at least for Frobenius manifolds of genus zero and one the Virasoro constraints 
hold true. That is WDVV equations can be reduced to that of KdV type. 
This result is in complete accord with arguments provided earlier in this section 
supporting the claim about universality of KdV type equations. This univer- 
sality has its origin in the properties of the universal Teichmiiller and moduli 
space. This universality has been studied group-theoretically by Schneps and 
Lochak [99] who studied action of the Grotendieck- Teichmiiller group GT on 
Ptolemy- Teichmiiller groupoid. Recent lecture notes by Bakalov and Kirillov 
[121 ] contain some additional helpful information. The authors also discuss 
connections between the Catalan numbers, modular functor, modular tensor 
category, Teichmiiller tower, Knizhnick-Zamolodchikov equations, etc. and con- 
tain many references on latest important related works. The latest paper by 
Schneps and Nakamura [122] should, perhaps, be added to the list references. 
It is our hope, that mentioning of all these beautiful mathematical results in 
this paper may encourage some applications of these results in areas other than 
mathematics and theoretical particle physics. 

Note added in proof.While this paper was under refereing, several impor- 
tant recent publications came to our attention. In particular, Ref.[123] provides 
a very efficient introduction into Schubert polynomials, Shubert varieties and 
related topics. Refs. [124,125] written and edited by Gian Carlo Rota provide 
an indispensable supplement to book by Stanley [42] . These referenses allow to 
keep things in historical perspective and, hopefully, should serve as an inspira- 
tion for many additional potential applications. In recent paper by Sozen and 
Bonahon [126] it is shown that the Weil Petersson symplcctic form ujw-p,e.g. 
see Eq.(1.4), coincides with Thurston intersection form t for geodesic lamina- 
tions ( more details on this subject can be found from yet unpublished book 
by Francis Bonahon "Closed Curves on Surfaces" available on line ).This allows 
to establish connections between foliations/laminations and non commutative 
geometry as actually had been demonstrated by Connes [127] some time ago. In 
addition to monograph [27] by Manin, recently Manin and Zograf had obtained 
explicit asymptotics for Weil-Petersson volumes of the moduli spaces of punc- 
tured Riemann surfaces in the limit of fixed genus and number of punctures 
groving to infinity [128] while Grushevsky have obtained the asymtotics for the 



48 



fixed number of punctures with genus groving to infinity [129]. In Ref. [130] 
Guha made some progress in showing new connections between the diffeomor- 
phisms of the circle and the exactly integrable systems of the Korteweg-deVries 
type by connecting them with differential Galois theory. Finally, direct links 
between the results of our Ref. [4] and those discussed in this paper could be 
also developed with help of random polynomials [131,132]. 
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